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Abstract 

We  present  a  queueing-theoretic  framework  to  evaluate  the  performance  of  large-scale,  query- 
based  wireless  sensor  networks  whose  nodes  detect  and  advertise  significant  events  that  are  useful 
for  only  a  limited  time;  queries  generated  by  the  nodes  of  the  network  are  also  time-limited.  The 
main  performance  parameter  is  the  steady  state  proportion  of  generated  queries  that  fail  to  be 
answered  on  time.  Using  an  infinite  transmission  range  model,  we  first  provide  an  approximation 
for  this  parameter  that  is  insensitive  to  the  size  of  the  network.  Subsequently,  we  approximate 
the  proportion  of  failed  queries  when  the  transmission  range  is  limited  and  show  that  this 
proportion  converges  to  its  infinite  range  counterpart  as  the  sensor  transmission  range  tends  to 
infinity.  The  analytical  approximations  are  shown  to  be  remarkably  accurate  when  compared 
with  benchmark  values  obtained  using  a  commercial  network  simulator. 


1  Introduction 

This  paper  evaluates  the  performance  of  large-scale,  query-based  wireless  sensor  networks 
(WSNs)  whose  sensors  detect  and  advertise  significant  events  that  are  useful  for  only  a  limited 
time  before  they  expire  (e.g.,  detecting  hazardous  biological  agents,  military  surveillance,  environ¬ 
mental  monitoring,  etc.).  Event  expiration  times  are  established  to  ensure  that  sensor  nodes  have 
the  most  up-to-date  information  to  share  with  other  nodes  in  the  network.  Query-based  WSNs 
derive  their  name  from  the  fact  that  communication  between  nodes  is  either  event-  or  query-driven. 
That  is,  either  the  witnessing  of  an  event  (e.g.,  a  sudden  increase  in  temperature),  or  the  generation 
of  a  query  (e.g.,  a  request  for  the  temperature  reading  at  a  distant  region  of  the  network)  triggers 
communication  between  nodes  which  must  act  as  routers  for  other  nodes’  packets  due  to  a  limited 
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sensor  transmission  range.  A  query,  which  itself  has  a  limited  lifetime,  traverses  the  network  accord¬ 
ing  to  a  two-dimensional  random  walk  until  it  either  locates  the  desired  information  or  expires.  For 
this  type  of  network,  the  total  proportion  of  generated  queries  that  are  not  answered  prior  to  their 
expiration  time  is  a  critical  performance  parameter.  We  devise  simple  analytical  approximations  for 
this  quantity,  along  with  other  network  quality-of-service  measures,  within  a  queueing  framework. 
Our  analytical  approach  is  unique  in  that  it  explicitly  accounts  for  the  realism  of  limited  event  and 
query  lifetimes  which  are  generally  distributed. 

Large-scale  wireless  sensor  networks  are  emerging  in  such  diverse  applications  as  ecological  and 
environmental  monitoring,  structural  health  monitoring  of  aging  infrastructure,  industrial  process 
control  and  military  surveillance.  The  ever-increasing  interest  in  WSNs  stems  from  their  ability 
to  sense  and  convey  critical  information  about  objects,  their  surroundings,  and  the  interaction 
between  the  two  autonomously.  Large-scale  WSNs  are  composed  of  thousands,  or  hundreds  of 
thousands,  of  low-cost  sensing  devices,  typically  linked  via  a  wireless  channel,  that  cooperate  to 
perform  specific  network  tasks  in  a  distributed  manner.  Due  to  their  small  physical  dimensions, 
the  sensing  nodes  have  very  limited  energy  reserves,  local  memory,  and  computational  capabilities. 
Moreover,  to  conserve  power  and  alleviate  contention  for  access  to  the  transmission  medium,  each 
node’s  transmission  range  may  be  limited  to  that  required  to  ensure  a  connected  network. 

Wireless  sensor  networks  have  been  analyzed  from  a  variety  of  perspectives  including  design 
considerations,  routing  protocols,  and  resource  management  strategies,  to  name  only  a  few.  Some 
useful  survey  papers  related  to  WSN  sensing  tasks,  applications,  design  issues,  and  communications 
architectures  include  Akyildiz  et  al.  [5],  Yick  et  al.  [37],  and  Dietrich  and  Dressier  [17].  Owing 
to  the  fact  that  sensor  nodes  are  energy-constrained,  defining  WSN  lifetime  and  operating  policies 
has  emerged  as  a  critical  issue.  Dietrich  and  Dressier  [17]  surveyed  many  definitions  of  WSN 
lifetime  including  the  number  of  “alive”  nodes,  network  coverage,  network  connectivity,  and  quality- 
of-service  considerations  (e.g.,  event  detection  rates).  Other  authors  (cf.  Anastasi  et  al.  [8]) 
have  classified  energy  conservation  approaches  (e.g.,  sensor  sleep/wake  protocols,  data  acquisition 
schemes,  mobile  sink-based  approaches,  etc.).  The  lifetime  of  WSNs  has  been  further  discussed  in 
[7,  28,  29,  38]. 

Routing  protocols  constitute  the  largest  area  of  research  related  to  the  performance  of  wireless 
sensor  networks.  A  variety  of  techniques  are  reviewed  in  a  cogent  survey  by  Al-Karaki  and  Karnal 
[6].  Most  protocols  aim  to  minimize  the  energy  expended  by  the  network  while  satisfying  quality- 
of-service  guarantees.  Routing  protocols  are  broadly  labeled  as  flat-based,  hierarchical  (see  [30]), 
and  location-based.  Flat-based  (or  data-centric)  protocols  assume  all  sensor  nodes  have  equal 
capabilities  and  similar  roles  whereas  hierarchical  protocols  assign  different  roles  to  the  nodes. 
Location-based  protocols  use  sensor  node  position  information  to  make  routing  decisions.  The 
model  we  analyze  here  falls  into  the  category  of  flat  routing  and,  more  specifically,  query-based  flat 
routing.  Classical  data-centric  approaches  include  flooding  and  gossiping  (see  [20])  which  are  known 
to  be  energy-  and  bandwidth- inefficient.  Alternatively,  rumor-routing  protocols  (see  [4,  14,  10,  18, 
31,  34,  36])  can  be  used.  Rumor  routing  uses  packets  with  relatively  long  lifetimes  called  agents. 
When  a  node  detects  an  event,  it  adds  information  pertaining  to  the  event  in  a  local  event  table 
and  immediately  creates  a  time-limited  agent  that  “advertises”  the  local  information  to  distant 
nodes  via  subsequent  packet  transmissions.  Consequently,  when  any  node  of  the  network  generates 
a  query,  a  node  with  the  information  stored  in  its  local  event  table  can  respond,  if  it  receives  the 
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query.  This  approach  obviates  the  need  for  flooding,  thereby  reducing  energy  expenditure.  Rumor 
routing  is  effective  (relative  to  flooding)  when  the  arrival  rate  of  events  is  relatively  small  but 
generally  requires  significant  overhead.  Specifically,  witnessed  events  are  assigned  a  time-to-live 
(TTL)  counter,  or  resource  replication  level ,  that  is  tracked  while  query  lifetimes  must  also  be 
tracked. 

The  resource  replication  level  is  the  number  of  times  a  witnessed  event  is  replicated  in  the 
network,  and  studies  related  to  this  parameter  are  relatively  sparse.  Bellavista  et  al.  [11]  developed 
a  simulation  model  (REDMAN)  to  explore  resource  replication  levels  and  related  network  settings. 
Krishnamachari  and  Ahn  [22]  derived  cost  expressions  as  a  function  of  the  resource  replication  level 
for  unstructured  networks  in  which  the  source  node  is  unknown.  They  used  expanding  ring  queries 
to  search  for  the  information  and  formulated  a  nonlinear  programming  (NLP)  model  to  determine 
the  optimal  number  of  resource  replicates,  subject  to  a  network  storage  capacity  constraint.  Ahn 
and  Krishnamachari  [3]  extended  the  results  of  [22]  to  structured  networks  in  d— dimensional  space, 
and  in  [1],  studied  structured  and  unstructured  two-dimensional  grid  and  random  topology  net¬ 
works.  Ahn  and  Krishnamachari  [2]  presented  a  model  to  obtain  the  optimal  resource  replication 
level  that  minimizes  the  total  expected  cost  of  replication  and  searching,  subject  to  a  storage  ca¬ 
pacity  constraint.  An  algorithm  for  dissemination  and  retrieval  of  information  that  ensures  an 
even  geographical  distribution  of  the  informed  nodes  is  proposed  for  unstructured  wireless  ad-hoc 
networks  by  Miranda  et  al.  [26].  Most  relevant  to  our  work  here,  Mann  et  al.  [25]  used  a  queueing 
framework  to  obtain  the  optimal  replication  level  that  minimizes  a  proxy  for  energy  expenditure, 
subject  to  a  performance  guarantee  on  the  steady  state  proportion  of  failed  queries.  Their  approach 
is  unique  in  that  it  considers  time-limited  event  agents  and  queries  but  is  limited  to  memoryless 
lifetimes.  Bisnik  and  Abouzeid  [13]  used  a  queueing  network  model  to  analyze  random  access, 
multi-hop  wireless  networks  and  derived  the  average  end-to-end  delay.  Niyato  and  Hossain  [27] 
developed  a  queueing  model  to  investigate  the  performance  of  different  sleep  and  wake-up  strate¬ 
gies.  Chiasserini  et  al.  [15]  proposed  a  fluid  queueing  model  that  accounts  for  energy  consumption, 
active/sleep  dynamics,  and  traffic  routing.  Ata  [9]  considered  the  problem  of  dynamically  choosing 
the  transmission  rate  in  a  general  wireless  communications  network  such  that  the  average  energy 
consumption  per  time  unit  is  minimized,  subject  to  a  quality-of-service  constraint.  In  that  work, 
the  transmission  queue  was  modeled  as  a  finite-buffer,  M/M/1  system.  With  the  exception  of 
Mann  et  al.  [25],  none  of  the  analytical  models  described  herein  account  explicitly  for  limited 
packet  lifetimes. 

This  paper  provides  a  queueing-theoretic  framework  for  evaluating  the  steady  state  proportion 
of  query  failures  (i.e.,  the  proportion  of  generated  queries  that  fail  to  be  answered  on  time)  in 
a  large-scale  WSN  with  time-critical  data.  While  the  network  model  itself  is  similar  to  the  one 
described  in  [25],  it  has  several  important  distinguishing  attributes.  Moreover,  the  analysis  ap¬ 
proach  we  present  here  is  quite  different.  In  fact,  the  model  of  Mann  et  al.  [25]  is  strictly  limited 
to  exponentially  distributed  inter-event  times  and  query  lifetimes,  and  it  essentially  ignores  the 
network’s  topology  by  assuming  that  sensors  have  an  infinite  transmission  range.  We  generalize 
their  model  by  considering  general  event  and  query  lifetimes  and  a  finite  transmission  range.  De¬ 
rived  are  valid  approximations  that  explicitly  account  for  (1)  time-limited  events  and  queries,  (2) 
the  limited  transmission  range  of  sensor  nodes,  and  (3)  generally-distributed  resource  and  query 
lifetimes.  The  first  approximation  (derived  from  an  infinite  transmission  range  model)  is  shown 
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to  be  insensitive  to  the  network’s  size.  The  finite  transmission  range  approximation  is  shown  to 
be  asymptotically  valid,  and  extensive  numerical  comparisons  with  simulated  networks  verify  the 
exceptional  accuracy  of  the  approximations,  even  when  exponential  assumptions  are  violated. 

The  remainder  of  the  paper  is  organized  as  follows.  Section  2  provides  a  description  of  the 
network  model,  the  queueing  models  of  sensor  node  elements,  and  the  most  relevant  attributes. 
In  Section  3,  we  derive  the  steady  state  proportion  of  query  failures  assuming  an  unlimited  sensor 
transmission  range,  while  Section  4  presents  an  approximation  that  explicitly  accounts  for  the 
limited  transmission  range  of  sensors.  Section  5  presents  extensive  numerical  results  that  validate 
our  analytical  approximations. 

2  Model  Description 

Consider  a  multi-hop  wireless  sensor  network  (WSN)  represented  by  an  undirected  graph  Q  = 
(AC  A)  where  N  =  {1,2, .. . ,  N}  is  the  node  set  (or  set  of  vertices),  N  is  the  number  of  sensor 
nodes  in  the  network,  and  A  is  the  arc  set  of  the  sensor  network.  An  arc  {%.  j)  is  an  element  of 
A  if  and  only  if  nodes  i  and  j  are  within  transmission  range  of  each  other.  Once  deployed,  the 
sensor  nodes  are  spatially  stationary  (i.e.,  they  are  not  mobile).  In  this  research,  we  consider  only 
networks  with  sensor  nodes  deployed  in  a  sensor  field  R,  a  subset  of  Euclidean  2-space.  The  nodes 
are  assumed  to  be  spatially  uniformly  distributed  in  R.  The  node  density  of  the  network,  ip,  is  the 
average  number  of  nodes  per  unit  area  (in  nodes/m2)  given  by  ip  =  N/A  where  A  is  the  area  of 
sensor  field  R.  For  each  i  £  A/",  denote  by  ccj  the  position  of  sensor  node  i  in  Euclidean  2-space. 
Then  for  j  £  A f ,  j  A  h  the  Euclidean  distance  between  and  Xj  is  p(i,j)  =  \\xi  ~  xj\\  where  ||  •  || 
denotes  the  Euclidean  norm.  Assuming  each  sensor  node  has  a  transmission  range  r  (in  meters), 
the  degree  of  node  i  £  M  is  the  number  of  nodes  within  transmission  range  of  i  given  by 

di(r)=  ^2  1(p(*.j)^r)> 

jeN\{i} 

where  l(x)  is  an  indicator  function  that  assumes  the  value  1  if  condition  x  holds  and  0  other¬ 
wise.  Obviously,  di{r)  depends  on  the  deployment  of  nodes  in  R,  the  network  topology,  and  the 
transmission  range  of  individual  sensor  nodes.  Finally,  the  average  degree  of  the  network  is 

1  N 

d(r )  =  ^2  di(r)- 

1=1 

A  node  i  £  AT  for  which  di(r )  =  0  is  said  to  be  isolated.  Isolated  nodes  are  essentially  useless 
to  the  WSN  since  they  cannot  exchange  information  with  other  nodes.  The  WSN  is  said  to  be 
disconnected  if  there  exists  at  least  one  isolated  node  in  J\f  but  is  completely  connected  if  there  exists 
at  least  one  path  between  nodes  i  and  j  for  every  i,  j  £  A f.  Obviously,  it  is  undesirable  for  the 
network  to  be  disconnected,  particularly  when  the  information  relayed  by  nodes  is  time  sensitive. 
When  the  nodes  are  uniformly  distributed  in  R  with  homogeneous  node  density  ip,  the  minimum 
transmission  range  needed  to  ensure  the  network  is  connected  with  probability  p  is  (see  Theorem 
1  of  [12]) 
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The  lower  bound  in  (1)  can  be  used,  for  example,  to  create  discrete-event  simulation  models  of 
wireless  sensor  networks  that  ensure  connectivity  with  high  probability. 

Next,  we  describe  individual  sensor  nodes  in  greater  detail.  (This  discussion  is  similar  to  that 
of  Mann  [25].)  It  is  assumed  that  sensor  nodes  are  identical,  i.e. ,  they  have  identical  resource 
requirements,  physical  limitations,  and  performance  limitations.  They  are  also  similar  with  respect 
to  their  information  requirements  and  the  rates  at  which  they  observe  and  report  relevant  phenom¬ 
ena.  Each  sensor  node  is  equipped  with  processing,  transmitting,  and  sensing  capabilities  as  well 
as  a  limited  power  supply  (in  the  form  of  an  on-board  battery)  that  is  generally  difficult,  if  not 
impossible,  to  replace. 

In  hybrid  push-pull  wireless  sensor  networks,  sensor  nodes  serve  as  both  producers  and  con¬ 
sumers  of  network  resources.  A  node  produces  a  resource  when  (1)  it  monitors  the  environment 
and  gathers  data  on  the  occurrence  of  pertinent  events;  or  (2)  it  offers  a  particular  service  to  the 
network.  In  addition  to  data  gathering,  nodes  are  also  required  to  execute  specific  applications 
in  support  of  the  network’s  goals.  When  a  node  requires  access  to  a  resource  that  is  not  avail¬ 
able  locally,  the  node  is  forced  to  traverse  the  network  to  locate  the  necessary  information  and/or 
services. 

When  a  node  witnesses  a  relevant  phenomenon,  or  offers  a  particular  service  to  the  network, 
it  broadcasts  this  information  to  a  subset  of  the  network  by  means  of  an  event  agent  -  a  packet 
that  describes  the  resource  available,  the  location  of  the  resource  (or,  alternatively,  the  data  itself), 
and  the  duration  of  time  the  resource  is  available  or  valid.  In  this  research,  we  assume  that  agents 
are  transmitted  from  node-to-node  via  a  random  walk  until  either  the  witnessed  event  expires 
(i.e.,  it  reaches  its  deadline),  or  it  exhausts  its  time-to-live  (TTL)  counter  -  a  positive  integer 
representing  the  maximum  number  of  times  the  resource  may  be  replicated  in  the  network.  It  is 
worth  mentioning  that  a  variety  of  routing  protocols  can  be  assumed  (cf.  [24]),  but  the  results 
herein  assume  transmission  to  a  randomly  selected  neighbor.  Each  sensor  node  is  equipped  with  an 
on-board  event  table.  Whenever  an  event  agent  is  received,  or  an  event  is  witnessed  by  the  node, 
the  contents  of  the  event  agent  are  added  to  the  event  table,  and  the  node  is  labeled  as  informed 
as  long  as  the  event  agent  has  not  expired. 

In  addition  to  witnessing  and  forwarding  events,  nodes  generate  queries  to  request  data  or 
resources  from  the  network.  A  query  contains  at  least  three  pieces  of  information:  the  identifier 
and/or  location  of  the  node  originating  the  request,  the  resource  sought,  and  the  maximum  amount 
of  time  the  query  is  permitted  to  traverse  the  network  in  search  of  an  informed  node.  Only  informed 
nodes  are  capable  of  answering  the  queries  of  uninformed  nodes.  Similar  to  event  agents,  queries 
are  forwarded  from  node-to-node  via  a  random  walk.  If  a  query  is  received  by  an  informed  node, 
the  query  is  terminated  and  the  informed  node  generates  a  response  that  is  returned  to  the  query 
origin  node  via  shortest-path  routing.  The  response  packet  contains  the  information  stored  in  the 
informed  node’s  event  table  and,  if  available,  the  desired  data.  If  a  query  cannot  locate  an  informed 
node  before  its  designated  expiration  time,  the  query  fails. 

Our  main  objective  is  to  assess  a  critical  quality-of-service  measure  for  query-based  WSNs, 
namely  the  long-run  proportion  of  queries  that  fail  to  be  answered  on  time.  To  this  end,  we  create 
a  queueing  network  model  that  leads  to  simple  analytical  expressions  and  accommodates  easy 
computational  implementation.  Next  we  describe  infinite-  and  single-server  queueing  models  of  the 
sensor  node’s  event  table  and  transmitter,  respectively. 
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2.1  Queueing  Models  of  Node  Elements 

For  each  i  £  J\f,  events  are  assumed  to  arrive  according  to  a  Poisson  process  with  rate  A.  Each 
witnessed  event  is  time  sensitive,  i.e.,  it  is  useful  for  only  a  limited  time  before  it  expires.  Therefore, 
once  an  event  is  witnessed  by  a  node,  it  is  added  to  the  node’s  event  table  and  assigned  an  expiration 
time  or  lifetime ,  Z,  a  non-negative,  non-defective  random  variable.  Event  expiration  times  (across 
all  nodes)  are  independent  and  identically  distributed  (i.i.d.)  random  variables  with  common 
cumulative  distribution  function  (c.d.f.)  G(w)  =  P(Z  <  w),  w  >  0,  and  mean  E (Z)  =  1/5  <  oo. 
As  long  as  the  event  agent  has  not  expired  in  the  event  table,  the  node  is  informed  and  can  answer 
queries  arriving  from  other  nodes.  Because  event  agents  are  mutually  independent  and  do  not 
necessarily  expire  in  their  order  of  arrival,  the  event  table  can  be  modeled  as  an  M/G/oo  queueing 
system  whose  input  is  a  Poisson  process  with  aggregate  arrival  rate  A  and  whose  service  time 
is  generally  distributed  with  c.d.f.  G  (see  Figure  1).  The  event  arrival  rate  A  depends  on  many 

Locally-witnessed  events 

Event  advertisements 

Event  expirations 

Figure  1:  Graphical  depiction  of  a  sensor  node’s  event  table  as  an  M/G/oo  queue. 

factors,  not  the  least  of  which  is  the  transmission  range  r.  We  pause  here  to  remark  that,  in  general, 
the  superposition  of  locally-witnessed  events  and  externally-generated  event  advertisements  does 
not  necessarily  form  a  Poisson  process  since  the  latter  do  not  (in  general)  originate  from  a  Poisson 
stream.  Furthermore,  the  event  table  may  not  realistically  have  infinite  capacity.  Therefore,  the 
evolution  of  the  number  of  busy  servers  in  the  M/G/oo  model  must  be  viewed  as  an  approximation 
of  the  evolution  of  event  table  content.  However,  it  will  be  shown  in  Section  5  that  this  assumption 
is  not  overly  restrictive  and  that  the  proportion  of  query  failures  is  surprisingly  insensitive  to  the 
Poisson  assumption.  We  choose  the  M/G/oo  model  for  its  tractability  and  generality  with  respect 
to  event  lifetimes.  Specifically,  it  provides  a  simple  expression  for  the  steady  state  proportion  of 
time  an  arbitrary  node  in  the  WSN  is  uninformed  given  by 

■kq  =  P(E  =  0)  =  exp  ( — A/ (5)  , 

where  E  is  the  steady  state  number  of  events  in  the  event  table.  Once  a  node  witnesses  an  event, 
the  information  is  forwarded  until  its  TTL  counter  is  exhausted.  Henceforth,  we  denote  this  integer 
value  by  t  £  N. 

Each  sensor  node  contains  a  transmitter  along  with  an  (assumed)  infinite  buffer  for  storing 
data  packets  (queries,  event  agents,  or  responses).  When  an  event  agent  arrives  to  a  node,  and  its 
time-to-live  counter  has  not  been  exhausted,  the  agent  joins  the  transmission  queue  after  a  copy 
has  been  added  to  the  node’s  event  table.  Moreover,  when  a  node  receives  a  query,  either  the  query 
or  the  response  is  sent  to  the  transmission  queue,  depending  on  whether  the  node  is  informed  or 
uninformed.  In  either  case,  the  query  fails  if  the  time  elapsed  from  the  moment  of  its  inception 
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until  it  locates  an  informed  node  exceeds  its  expiration  time.  The  node’s  transmission  queue  is 
modeled  as  a  single-server  queueing  system  as  depicted  in  Figure  2. 

Locally-generated  queries 
External  queries 
Locally-witnessed  events 
Event  advertisements 

vV 

Event  and  query  expirations 

Figure  2:  Graphical  depiction  of  the  sensor  node’s  transmission  queue. 

Specifically,  we  assume  that  each  node’s  transmission  queue  operates  as  a  non-prioritized,  multi¬ 
class  M/M/1  queueing  system  with  a  first-come- first-served  (FCFS)  queueing  discipline.  Arrivals 
are  assumed  to  originate  from  a  Poisson  process  with  rate  Xq.  As  depicted  in  Figure  2,  the  aggregate 
arrival  process  is  comprised  of  locally-witnessed  events,  advertisements  from  other  nodes,  locally- 
generated  queries,  and  queries  arriving  from  other  nodes.  When  a  query  is  generated  or  received 
by  a  node,  it  joins  the  transmission  queue  only  if  the  node  is  uninformed.  The  service  time  is  the 
time  needed  to  transmit  a  query  packet  or  an  event  packet  (either  a  locally-witnessed  event  or  an 
advertisement  from  other  nodes).  Irrespective  of  the  packet  type,  we  assume  the  transmission  time 
r  is  an  exponential  random  variable  with  parameter  fj,,  c.d.f.  F(x)  =  P(r  <  x)  =  1  —  exp(— fix),  and 
finite  mean  E(r)  =  1/fj,.  The  transmission  queue  is  stable  if  and  only  if  ji  >  Xq.  This  condition  is 
usually  met  in  practice  since  transmission  rates  are  generally  very  high.  It  is  important  to  note  that 
the  total  arrival  rate  of  traffic  to  the  transmission  queue  serves  as  a  proxy  for  energy  expenditure 
at  a  node  since  transmitting  is  the  primary  energy  consuming  activity  (cf.  [32]). 

2.2  Network  Performance  Parameters 

The  primary  concern  of  this  research  is  the  assessment  of  the  steady  state  probability  that 
a  generated  query  fails  to  be  answered  on  time.  We  refer  to  this  performance  parameter  as  the 
proportion  of  query  failures.  A  query  is  said  to  fail  if  it  expires  awaiting  transmission  or  while  being 
transmitted.  Our  main  aim  is  to  provide  easy-to-use  analytical  expressions  for  this  parameter  that 
allow  us  to  circumvent  costly,  time-consuming  simulation  experiments  for  large-scale  networks.  To 
this  end,  let  T  be  a  non-negative  random  variable  denoting  the  total  time  needed  for  a  query  to 
locate  an  informed  node  as  measured  from  the  time  the  query  is  generated  at  a  node  n  £  M .  This 
random  time  depends  on  the  status  of  node  n  at  the  time  of  creation.  Define  the  indicator  variable 

1,  if  node  n  is  informed, 

0,  if  node  n  is  uninformed. 

The  c.d.f.  of  [T\In  =  0]  is  denoted  by  B(t)  =  P(T  <  t\In  =  0),  t  >  0  for  any  n  £  N .  At  the  time  of 
generation,  the  query  is  assigned  a  lifetime,  X ,  so  that  if  T  exceeds  X ,  the  query  does  not  locate 
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the  desired  information  before  expiring,  and  it  fails.  The  c.d.f.  of  X  is  H(x)  =  P(A  <  x ),  x  >  0, 
and  its  mean  is  E(X)  =  1//3  <  oo.  Proposition  1  characterizes  the  primary  performance  parameter. 

Proposition  1  The  unconditional  proportion  of  query  failures  is 

roo 

A  =  P(T  >  X)  =  tto  /  [1- B(x)]dH(x).  (2) 

Jo 

Proposition  1  can  be  proved  using  a  simple  conditioning  argument.  The  expression  for  the 
proportion  of  query  failures  is  straightforward  except  that  the  distribution  function  B  is  difficult 
to  characterize  in  all  but  a  few  cases.  The  next  section  considers  one  such  case  when  the  sensors 
all  use  an  infinite  transmission  range. 

3  Unlimited  Sensor  Transmission  Range 

In  this  section,  we  provide  an  approximation  for  A  when  r  =  oo,  i.e.,  when  each  node’s  trans¬ 
mission  range  is  large  enough  to  ensure  that  any  other  node  in  the  network  can  be  reached  with 
a  single  hop.  This  simplified  model  provides  an  approximation  for  the  proportion  of  failed  queries 
because  realistic  sensor  nodes  have  limited  transmission  capabilities  and  rely  on  their  neighbors  to 
propagate  witnessed  events,  queries,  and  responses  to  queries.  Although  several  assumptions  are 
employed,  the  primary  purpose  of  this  model  is  to  provide  a  framework  for  a  more  realistic  limited 
transmission  range  model. 

3.1  Approximating  Network  Traffic 

Here  we  establish  approximations  for  event  agent  and  query  arrival  rates  at  an  arbitrary  node 
of  the  network,  assuming  nodes  are  identical.  The  first  result  bounds  the  aggregate  event  arrival 
rate  to  the  sensor  node’s  event  table.  This  bound  sets  the  stage  for  an  approximation  of  the  steady 
state  proportion  of  time  that  any  node  is  uninformed. 

Proposition  2  Assume  events  arrive  locally  to  each  n  £  J\f  according  to  a  Poisson  process  with 
rate  A.  Then  the  total  event  arrival  rate  A  to  an  arbitrary  n  £  J\f  is  bounded  above  by  A  (1  +  t) 
where  i  is  the  time-to-live  counter. 

Proof.  The  aggregate  event  arrival  rate  A  consists  of  the  Poisson  rate  of  locally-witnessed 
events,  and  the  aggregate  rate  of  witnessed  events  arriving  from  the  other  N  —  1  nodes  in  the 
network.  Therefore,  A  =  A  +  A^  where  Ax  denotes  the  rate  of  external  event  arrivals.  An  event 
agent  can  be  forwarded  to,  at  most,  l  nodes.  Since  r  =  oo,  each  node  can  be  reached  with  a  single 
hop;  hence,  each  event  advertisement  is  equally  likely  to  be  received  by  one  of  the  other  N  —  1 
nodes.  That  is,  a  particular  node  receives  one  of  the  (potential)  i  advertisements  with  probability 
£/(N  —  1),  and  since  N  —  1  other  nodes  transmit  event  agents, 

A*<A(IV  — 1)  (^-)  =A1 


Therefore,  A  <  A  +  \l  =  A  (1  +  l). 


Next,  we  provide  a  lower  bound  for  the  steady  state  proportion  of  time  a  node  is  uninformed. 
An  event  agent  is  assigned  a  lifetime  Z  once  it  enters  the  event  table.  The  random  variable  Z 
has  c.d.f.  G  and  finite  mean  E (Z)  =  1/5.  As  noted  in  Section  2,  the  event  table  is  approximated 
by  an  Al/G/oo  queue  with  (Poisson)  arrival  rate  A  due  to  the  insensitivity  of  this  model  to  the 
event  lifetime  distribution.  Using  well-known  results  for  the  steady  state  behavior  of  the  M/G/oo 
system,  the  limiting  proportion  of  time  a  node  is  uninformed  (i.e.,  the  limiting  probability  that  the 
event  table  is  empty),  is 

7To  =  exp(— A/<5),  0  <  6  <  oo. 

By  Proposition  2,  the  event  arrival  rate  to  a  node  is  bounded  above  by  A  (1  +  £).  Therefore, 


7Tq  >  exp 


a(i  +  ey 
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(3) 


Similarly,  the  event  arrival  rate  to  the  node’s  transmission  queue ,  Ae,  is  also  bounded  above. 
Specifically,  Ae  <  X£  since  a  locally-witnessed  event  is  certainly  added  to  the  transmission  queue. 
On  the  other  hand,  when  a  witnessed  event  is  advertised  through  the  network,  since  each  of  the 
other  N  —  1  nodes  is  equally  likely  to  receive  the  advertisement,  a  single  node  receives  the  first 
transmission  and  adds  the  event  agent  to  its  transmission  queue  with  probability  1  /{N  —  1),  the 
second  transmission  with  the  same  probability,  and  so  forth.  The  event  agent  is  forwarded  at 
most  £  times,  but  the  node  which  receives  it  at  the  fth  transmission  does  not  add  the  agent  to  its 
transmission  queue  since  the  agent’s  time-to-live  counter  will  have  expired.  Therefore, 
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Ae<A  +  A(AT-l)^^-I  =  AU  (4) 

While  the  bounds  of  (3)  and  (4)  are  legitimate,  they  may  not  be  tight  since  they  do  not  explicitly 
account  for  the  expiration  of  event  agents  waiting  in  the  transmission  queue,  or  those  that  expire 
during  transmission.  Proposition  3  provides  an  improved  bound  for  Ae  (by  considering  the  effect  of 
event  expirations)  that  leads  to  an  improved  approximation  for  7Tq.  In  what  follows,  let  aj  denote 
the  probability  that  an  event  agent  expires  at  the  jth  visited  node,  j  =  1,2,...,  where  the  first 
visited  node  is  the  event  witnessing  node.  For  simplicity,  define  the  expiration  probability  at  the 
first  node  by  a  =  ot\.  Assuming  the  event  lifetime  distribution  function  G  has  an  increasing  failure 
rate  (IFR),  then  0  <  a  <  «2  <  0:3  <  ■  ■  ■ . 


Proposition  3  Suppose  G  is  an  IFR  distribution  function  so  that  0  <  a  <  «2  <  «3  <  •  •  • .  Then 
for  a  fixed  time-to-live  counter  £, 


Ae  <  A 


'1  -  (1  -  ay 

a 


<  XL 


Proof.  Since  each  of  the  N  —  1  nodes  is  equally  likely  to  receive  an  advertised  event  agent,  an 
individual  node  receives  the  kth  transmission  with  probability 

k 

k  =  l,2,...,£-l. 

3= 1 
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That  is,  the  event  agent  is  forwarded  at  most  t  times;  however,  the  last  node  that  receives  the  agent 
does  not  add  it  to  its  transmission  queue  since  the  agent’s  time-to-live  counter  will  have  expired. 
Therefore,  the  approximate  total  rate  of  event  arrivals  to  a  node’s  transmission  queue  is  given  by 

t—\  k 

A«  =  A  +  AtJV-D^^h-nil-a,) 

k= 1  j= 1 

e- 1 

<  A  +  A-^(l-a)fc 

k=\ 

i-l 

=  A£(l-a)‘. 

k= 0 

Since  a  €  (0, 1),  we  apply  the  well-known  identity  for  geometric  summations 


E 

fc=0 


1  —  x 


m+1 


X  = 


1  —  X 


x  €  [0, 1), 


to  obtain  the  final  result, 


AP  A 


1  -  (1  -a)1 


a 


<M,  a  €(0,1), 


where  the  second  inequality  follows  from  a  €  (0, 1). 


For  the  results  that  follow,  we  use  the  approximation, 


A,  «  A 


1  -  (1  -a) 


a 


to  improve  the  approximation  of  7Tq.  Similarly,  it  can  be  shown  that  the  event  arrival  rate  to  the 
event  table  is  approximated  by 


A 


A  =  A 


'i  -  (i  -  ay+i' 

a 


(5) 


Therefore,  when  r 
given  by 


oo,  the  approximate  steady  state  proportion  of  time  a  node  is  uninformed  is 


7T0 


exp 


a  /i  -  (i  -ay+i 

8  \  a 


(6) 


Next,  we  examine  the  total  traffic  experienced  at  the  transmission  queue.  Let  Xq  be  the  total 
arrival  rate  of  events  and  queries  to  a  sensor  node’s  transmission  queue.  Each  node  generates  local 
queries  according  to  a  Poisson  process  with  rate  7.  When  a  query  is  generated  locally,  or  received 
from  another  node,  it  is  added  to  the  transmission  queue  only  if  the  subject  node  is  uninformed. 
The  arrival  rate  of  locally-generated  queries  to  the  transmission  queue,  qi ,  is  qi  =  vroy.  Let  qx 
denote  the  rate  at  which  external  queries  arrive  at  a  node.  In  steady  state,  the  query  visits  an 
informed  node  with  probability  1  —  ttq  and  an  uninformed  node  with  probability  7To,  independent 
of  the  number  of  hops  prior  to  the  current  visit.  Consequently,  the  number  of  hops  before  a  query 


10 


locates  an  informed  node  is  a  geometric  random  variable  with  success  probability  1  —  ttq,  and  mean 
1/(1  —  7To).  Since  any  node  is  equally  likely  to  receive  a  query,  the  probability  of  receiving  a  query 
generated  by  any  other  node  is  1/((1  —  ttq)(N  —  1)).  Because  all  other  N  —  1  nodes  generate  queries 
identically,  qx  is  approximately 


Qx  =  7^0  (N  -  1) 


1 


L  (i —  7to)(n  —  i) 


7^0 

1  -  7T0' 


Consequently,  we  approximate  the  total  arrival  rate  of  traffic  to  a  node’s  transmission  queue  as 


An  ~  A e  +  qi  +  qx  —  X 


1  -  (1  -af 


a 


+  7T07 


2  —  7Tq 
1  -  7T0 


(7) 

(8) 


By  (6)  and  (8),  we  see  that  -kq  and,  consequently,  \q  are  explicit  functions  of  a.  Therefore,  the 
approximation  of  \q  is  written  as 


Aq 


c(a )  =  A 


1  -  {1-a.y 


a 


+  ye 


-9(a) 


2  -  e~3(a) 

1  _  e-g(a)  ’ 


(9) 


where 


5(a) 


A  r  1  —  (1  -a)i+1' 
6  a 


We  are  now  prepared  to  provide  an  expression  for  a,  the  probability  that  an  event  agent  expires 
in  the  first,  transmission  queue.  The  result  is  approximate  since  the  input  to  the  transmission  queue 
is  assumed  to  be  the  superposition  of  independent  Poisson  arrival  streams.  The  equilibrium  random 
variable  Ze  associated  to  the  lifetime  Z  with  c.d.f.  G  and  mean  E (Z)  has  c.d.f. 


Ge(z)  =  P(Ze  <Z)  =  J^J~[  1  -  G(u)]du. 

We  make  use  of  the  equilibrium  distribution  of  the  event  agent  lifetime  in  the  following  proposition 
that  characterizes  a. 


Proposition  4  Assume  p,  >  c(a )  for  each  A  and  5  such  that  0  <  A  <  oo  and  0  <  5  <  oo.  Let  W  be 
the  total  time  spent  at  a  node ’s  transmission  queue  ( delay  plus  transmission  time )  by  an  arbitrary 
arrival  in  steady  state.  Then 


a  =  P (W  >  Ze) 


Ge(fr 


c(a ))  =  E 


e-\u~c(a)]Ze 


(10) 


where  Ge(p  —  c(a))  denotes  the  Laplace- Stieltjes  transform  (LST)  of  Ge  evaluated  at  p  —  c(a). 


Proof.  The  transmission  queue  is  modeled  as  an  M/M/1  queueing  system  with  mean  transmis¬ 
sion  time  1  / [i  and  aggregate  arrival  rate  c(a).  Let  Wn  be  the  total  time  spent  in  the  transmission 
queue  (i.e. ,  the  delay  time  plus  the  transmission  time)  by  the  nth  arrival  to  the  queue,  either  an 
event  agent  or  a  query.  It  is  well  known  that  if  //  >  c(a),  Wn  =4>  W  as  n  — »  oo  where  W  is 
exponentially  distributed  with  mean  l/(n  —  c(a ))  and  (=>)  is  convergence  in  distribution  (or  weak 
convergence).  Suppose  an  event  agent  arrives  at  time  t  so  that  Z  —  t  is  the  residual  lifetime  of  the 
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event.  Using  basic  results  from  renewal  theory,  Z  —  t  =$■  Ze  as  t  — >  oo.  Therefore,  by  conditioning 
on  Ze,  we  obtain 


a  =  P (w  >  Ze)  =  /  e-(M-c(a)Ud Ge(z)  =  Ge( [I  -  c(a))  =  E 


0-(/i-c(a))Ze 


As  an  illustration,  if  the  event  lifetime  is  exponentially  distributed  with  mean  1/5,  then  Ge(z)  = 
G(z)  for  all  z  >  0,  and  the  unique  probability  a  is  the  solution  to  the  fixed  point  problem 

5 

a  = - — ?• 

H  —  c(a)  +  o 

Recall  that  our  aim  is  to  approximate  A  of  equation  (2)  by  assuming  r  =  oo.  To  this  end,  let 
T  denote  the  time  to  locate  an  informed  node  when  r  =  oo  and  let 

~  r°° 

Aoo  =  P(T  >  X)  =  7T0  /  [1  -  B(x)]dH(x). 

Jo 

The  c.d.f.  of  T  is  a  function  of  both  Xq  (c(a))  and  ttq,  both  of  which  are  determined  by  a.  The 
next  section  shows  how  to  obtain  Aoo. 

3.2  Approximate  Query  Failure  Rate 

Queries,  which  can  be  generated  at  any  node  n  €  A/",  are  forwarded  via  a  random  walk  to 
one-hop  neighbors  until  either  an  informed  node  is  located  or  the  query  expires  while  awaiting 
transmission  in  some  node’s  transmission  queue  (or  while  being  transmitted).  Once  generated,  a 
query  is  assigned  a  lifetime  X  having  c.d.f.  H(x)  =  P(A  <  x),  x  >  0.  Let  us  assume  for  the 
moment  that  a  query  generated  at  an  uninformed  node  can  be  forwarded  indefinitely  (i.e.,  X  =  oo 
w.p.  1),  and  let  M  be  the  (integer)  number  of  hops  needed  to  first  locate  an  informed  node. 

Let  Tk  denote  the  time  spent  by  a  query  at  its  /cth  location.  That  is,  T\  denotes  the  time  spent 
at  the  original  query  node  (which  is  uninformed),  T2  is  the  time  spent  at  the  second  node,  which 
might  be  informed  or  uninformed,  and  so  forth.  To  simplify  notation,  let  Tu  =  [T \ In  =  0]  be  the 
elapsed  time  between  creation  of  the  query  at  an  uninformed  node  and  the  time  it  first  locates  an 
informed  node.  It  is  easy  to  see  that 

M 

tu  =  y^Tfc. 

k= 1 

Because  we  assume  r  =  00  and  identical  nodes,  a  query  visits  an  informed  node  with  probability 
1  —  ttq  and  an  uninformed  node  with  probability  7To,  independent  of  any  prior  visits,  in  steady  state. 
Thus,  M  is  a  geometric  random  variable  with  success  probability  1  —  ttq  and  mean  1/(1  —  7To),  i.e., 
Tu  is  a  geometric  sum  of  i.i.d.  exponential  random  variables. 

Lemma  1  Given  that  a  query  is  generated  at  an  uninformed  node,  the  time  to  locate  an  informed 
node  is  exponentially  distributed  with  parameter  (1  —  7ro)(/U  —  \q),  i.e., 

B(x)  =  P(T„  <  x)  =  1  -  exp  [-(1  -  7T0)(/U  -  Ag)x] ,  x  >  0,  (11) 

where  \q  =  c(a )  is  obtained  using  the  value  of  a  that  solves  the  fixed  point  equation  (10). 
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Using  Lemma  1,  we  next  provide  our  approximate  expression  for  the  steady  state  proportion  of 
query  failures  when  r  =  oo. 

Proposition  5  Assuming  Poisson  event  arrivals  and  query  generation,  the  proportion  of  query 
failures  in  an  infinite-range  WSN  is 


Aoo  =  P(T  >X)  =  7T0  H[(  1  -  7T0)(fi  -  A,)],  (12) 

where  H(s)  =  E  (e_s*Y)  is  the  LST  of  the  query  lifetime  distribution  function  H. 

Proof.  The  proof  follows  directly  by  conditioning  on  the  lifetime  X  and  utilizing  Lemma  1. 
Specifically, 


Aoo  =  P(T  >  X) 


P(T>  X\In  =  0,X  =  x)P{In 

roo 

7 r0  /  e-U-^oXM-A^d H(x) 

Jo 


tt0H[(1  -  n0)(n  -  Aq)]  . 


O)dH(x) 


Proposition  5  holds  for  all  query  lifetime  distributions  that  possess  a  LST.  However,  if  the  distribu¬ 
tion  function  H  is  heavy-tailed  and  does  not  possess  an  LST,  the  transform  approximation  method 
(TAM)  developed  by  [19],  or  its  modification  by  [35],  can  be  used  to  approximate  H. 

3.3  Computing  the  Proportion  of  Query  Failures 

Here,  we  describe  a  simple  fixed  point  algorithm  to  compute  the  probability  a  to  obtain  A^ 
using  equation  (12).  Let  \'q  -  and  a ^  be  the  approximated  values  of  7To,  Xq  and  a  at  the  kth 
iteration  of  the  algorithm,  respectively. 


Algorithm  to  Compute  a: 

Step  0:  Initialization  via  the  bounds  of  (3)  and  (4). 
k  :=  0; 

n(0k)  ■.=  exp  [— A(1  +  f)/h]; 


;=  A  ^  +  77Tofe)  (  - - — 


1  -  7T, 


(fc)' 

0 

(fc) 


:=  Ge 


Step  1:  Update  the  approximations, 
k  :=  k  +  1; 

A  / 1-  (!  - 


( k ) 

7Tq  :=  exp 


<+i\  n 


Af  :=  A 


1_(1_  „(*-!)) 


a(k  1) 

to 


a 


(fc— i) 


(U 

+  7^0 


2  —  7T, 


(fc) 

0 


1  —  7T, 


(fc) 

0  J 
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Step  2  :  Check  convergence  criterion. 

If  |  >  e,  return  to  Step  1; 

Else  a  :=  a^; 

Stop. 

Scalability  of  the  WSN  is  an  important  issue  as  realistic  networks  are  envisioned  to  have  thou¬ 
sands  or  even  hundreds  of  thousands  of  sensor  nodes.  The  infinite  range  approximations  of  this 
section  are  appealing  due  to  their  insensitivity  to  N.  For  large  N,  the  probability  that  a  given 
node  is  visited  more  than  once  by  a  specific  event  agent  or  query  is  negligible.  That  is,  when  a 
node  witnesses  an  event,  it  forwards  an  event  packet  to  at  most  t  distinct  nodes.  Similarly,  when  a 
query  is  generated,  the  present  model  assumes  it  visits  a  distinct  node  at  each  hop,  independently 
of  all  prior  hops.  However,  to  conserve  power,  realistic  WSN  sensors  use  a  limited  transmission 
range,  so  the  probability  of  revisiting  neighbors  can  be  significant  as  highlighted  by  Rodero-Merino 
et  al.  [33].  This  revisiting  effect  increases  the  proportion  of  time  nodes  are  uninformed,  the  time 
to  locate  an  informed  node,  and  consequently,  the  proportion  of  failed  queries.  In  the  next  section, 
we  account  for  a  limited  sensor  transmission  range. 

4  Accounting  for  Transmission  Range 

In  this  section,  we  present  an  approximation  for  the  steady  state  proportion  of  query  failures 
that  explicitly  accounts  for  the  limited  transmission  range  of  wireless  sensors.  Specifically,  we 
approximate  Ar,  the  steady  state  proportion  of  query  failures  when  the  sensor  nodes  have  trans¬ 
mission  range  of  r  (r  <  oo).  Additionally,  we  show  that  for  large  N,  the  approximation  converges 
appropriately  to  as  r  — >•  oo. 

4.1  Modeling  Query  Dynamics 

Here  we  consider  the  status  (and  movement)  of  an  individual  query  from  its  inception  until  it 
either  locates  an  informed  node  or  fails  due  to  expiration.  If  a  query  is  generated  at  an  informed 
node,  then  it  is  answered  immediately  and  never  forwarded;  therefore,  we  focus  on  the  case  when 
a  query  is  generated  at  an  uninformed  node  n  G  A f.  At  its  inception  the  query  is  instantaneously 
assigned  a  lifetime  (or  expiration  time)  X  with  c.d.f.  H(x)  and  mean  E(A)  =  1//3  (0  <  /3  <  oo). 
It  is  forwarded  to  a  randomly  selected  node  within  the  ?’-radius  of  the  current  node  until  either  an 
informed  node  is  located,  or  the  query  expires,  in  which  case  it  is  destroyed.  In  what  follows,  all 
random  quantities  are  conditioned  on  the  event  {X  =  ir};  therefore,  we  make  the  dependence  on  x 
explicit. 

For  each  integer  k  (k  >  0),  let  Qk  be  the  status  of  the  query  just  before  it  joins  the  transmission 
queue  of  the  fctli  visited  node.  (Note  that  Q o  is  the  query’s  status  at  the  moment  it  is  generated 
at  the  query  origin  node.)  The  query  can  be  in  one  of  three  mutually  exclusive  and  exhaustive 
states:  active  (state  0),  answered  (state  1),  or  expired  (state  2).  The  query  is  added  to  the  kth 
visited  node’s  transmission  queue  if  and  only  if  it  is  active.  For  each  k  >  0,  Qk  G  S  =  {0, 1,2} 
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where  Qk  =  0  means  the  query,  having  been  successfully  transmitted  k  times,  has  not  expired  but 
has  not  been  answered;  Qk  =  1  means  the  query,  having  been  successfully  transmitted  k  times 
is  answered  at  the  kth  visited  node  (i.e.,  the  kth  visited  node  is  informed);  and  Qk  =  2  means 
the  query  was  successfully  transmitted  k  —  1  times  but  expired  awaiting  its  kth  transmission  (or 
during  its  kth  transmission)  at  the  (k  —  l)st  visited  node.  (Note  that  P(Qo  =  2)  =  0.)  The 
process  Q  =  {Qk  :  k  >  0}  is  an  5'- valued  discrete-time  Markov  chain  (DTMC)  with  temporally- 
nonhomogeneous  one-step  transition  probability  matrix,  P (k,x),  given  by 

p00(k,x)  poi{k,x )  p02(k,x) 

P (k,x)  =  0  1  0  ,  /c  >  0,  x>0  (13) 

0  0  1 

where  for  each  i.  j  £  S, 

Pij(k,x )  =  P.T(Qfc+i  =  j\Qk  =  i),  k>t), 

denotes  the  probability  that  the  status  of  the  query  transitions  from  i  to  j  at  the  (k  +  l)st  step  and 
PX(T)  =  P(j4|X  =  x)  for  any  measurable  event  A.  Once  a  query  locates  an  informed  node,  it  is 
no  longer  forwarded  to  a  neighbor  node,  and  if  the  query  expires  awaiting  transmission  (or  during 
its  transmission),  it  is  destroyed;  therefore,  states  1  and  2  are  absorbing  states  of  the  DTMC.  Row 
0  of  P (k,x)  contains  the  critical  transition  probabilities.  In  particular,  po2(k,x)  is  the  probability 
that  the  query  fails  at  the  kth  visited  node,  given  it  was  active  just  before  being  added  to  the  kth 
visited  node’s  transmission  queue.  Likewise,  poo(k,x)  is  the  probability  the  query  remains  active 
just  before  being  added  to  the  (k  +  l)st  visited  node’s  transmission  queue,  given  it  was  active  just 
before  being  added  to  the  kth  node’s  transmission  queue.  Finally,  poi(k,x)  is  the  probability  that 
a  query  is  answered  at  the  (k  +  l)st  visited  node,  given  it  was  active  just  before  being  added  to  the 
/cth  visited  node’s  transmission  queue. 

Obviously,  the  DTMC  {Qk  '■  k  >  0}  is  reducible  with  one  transient  state  (state  0)  and  two 
closed  communicating  classes,  namely  C\  =  {1}  and  C2  =  {2};  therefore,  its  limiting  behavior  is 
fairly  easy  to  characterize.  Before  examining  the  limiting  behavior,  we  characterize  the  distribution 
of  {Qk  :  k  >  0}  at  a  particular  step  k.  Let  vk(x)  =  P x{Qk  =  j)  be  the  (unconditional)  probability 
that  the  query  is  in  state  j  £  S  just  before  joining  the  transmission  queue  of  the  kth  node,  and 
let  vk(x)  =  [vk(x)]j£s  be  a  (1  x  3)  row  vector  comprised  of  these  values.  Because  {Qk  :  k  >  0} 
possesses  a  time-nonhomogeneous  transition  probability  matrix,  the  vector  vk(x)  can  be  obtained 
recursively  (cf.  [21])  by 

vk+1(x)  =  vk(x)  P(/c,  x),  k  >  0, 

whose  solution  is  given  by 

k 

vk+1(x)  =  v°(x)  P(n,  x),  k>  0.  (14) 

71=0 

The  square  matrix  on  the  right-hand  side  of  (14)  is  the  (/c  +  l)-step  transition  probability  matrix 
of  {Qk  '■  k  >  0}.  The  transient  analysis  of  {Qk  ■  k  >  0}  will  facilitate  an  analysis  of  its  limiting 
behavior  which,  in  turn,  is  used  to  derive  an  expression  for  the  steady  state  probability  that  a  query 
fails  to  locate  an  informed  node  before  expiring. 
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To  this  end,  let  us  define  the  limiting  probability  vector 

k  k 

v{x)  =  lim  vk+1(x )  =  lim  v°(x)  TT  P  (n,x)  =  v°(x)  lim  ][  P  (n,x).  (15) 

k—> oo  k — ^oo  k—> oo 

n=0  n= 0 

Before  approximating  this  vector,  we  first  establish  the  existence  and  form  of  the  limit  in  the 
right-most  term  of  (15)  via  Theorem  1. 


Theorem  1  For  a  fixed  expiration  time  x  (x  >  0),  there  exist  real  numbers  a i(x)  and  0:2(2;)  such 
that 

k 

A(x)  =  lim  1  I  P (n,x)  = 

k-yoo  1  x 
n= 0 

where  a.\{x),  02(2;)  G  (0,1)  and  Y^i= \ai{x)  =  1- 


0  «i(x)  02(2;) 

0  1  0 
0  0  1 


Proof.  Using  induction,  it  can  be  shown  that  the  (k  +  l)-step  transition  probability  matrix  is 
given  by 


k 

n  p(n^) 


n= 0 


nto  Etc  K  (n”;„‘  Oj)  Eto  On  «.,) 

0  1  0 

0  0  1 


(16) 


where  an  =  poo(n,x),  bn  =  po\{n,x),  cn  =  po2{n,x),  n  >  0,  and  a_i  =  1.  First,  note  that  rows  1 
and  2  of  n!)=oP(n>x)  are  as  giyen  in  (16)  for  any  k  €  N;  hence,  we  need  only  concern  ourselves 
with  row  0.  Allowing  k  — >  00  on  both  sides  of  (16),  and  noting  that  0  <  an  <  1,  we  see  immediately 


that 


k 


n= 0 


k 

n— 1 

00 

n—  1 

a\{x) 

=  lim 

k^f-oo 

Eb- 

n°j 

=  ^0  +  bn 

aj  >  b0  >  0, 

(17) 

n= 0 

3= 0 

n=  1 

j= 0 

k 

n— 1 

00 

71—1 

o2(x) 

=  lim 

k—>  00 

n  aj 

=  co  +  cn 

aj  >  co  >  0. 

(18) 

n= 0 

j= 0 

n=  1 

i=o 

Since  each  row  of  A(x)  is  comprised  of  nonnegative  real  numbers,  and  the  row  sums  must  be 
unity  (cf.  [16]),  we  conclude  that  a\(x)  +  a.2(x)  =  1  which,  in  light  of  (17)  and  (18),  implies  that 
0  <  0:1(2:)  <  1  and  0  <  02  (x)  <  1.  ■ 


For  computational  purposes,  we  approximate  v(x)  by  truncating  the  infinite  product  of  (15)  at 
an  appropriate  integer  q.  Specifically,  for  a  sufficiently  large  q  £  N,  the  approximation  for  v(x)  is 
given  by 

q 

v(x)  «  vq+1(x)  =  v°(x)  J^[  P(n,x),  (19) 

n= 0 
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where  q  is  chosen  such  that  ||v9+1(a;)  —  vq(x)\\00  <  e  with  ||  •  Hoc  the  usual  oo-norm  and  e  a 
convergence  threshold.  Let  Ar  be  the  limiting  probability  of  query  failure  provided  each  sensor’s 
transmission  range  is  r  (r  <  oo)  and  let 

V2(x)  =  lim  v%(x)  =  lim  P x(Qk  =  2). 

k—>  oo  k—>  OO 

The  unconditional  proportion  of  query  failures  is  approximately 

roo 

A  r  =  /  V2(x)d  H(x).  (20) 

Jo 

Let  7To (r)  be  the  steady  state  proportion  of  time  an  arbitrary  node  is  uninformed  when  the 
transmission  range  is  r.  Note  that  v^{x)  depends  implicitly  on  r  through  7To(?’)  since  u°(x)  = 
(7To(?’),  1  —  vro(r),  0),  and  A(x)  depends  on  7To(r).  However,  we  suppress  this  dependence  on  r 
for  ease  of  notation.  To  compute  v(x)  (or  its  approximation  vq(x)  via  (19)),  we  now  provide  an 
expression  for  po2{k,x)  and,  subsequently,  expressions  for  poo(k,x)  and  po\(k,x). 


Lemma  2  For  a  fixed  expiration  time  x  (x  >  0),  the  transition  probability  po2(k,  x)  is 

[(yU  —  A  q)x]k 


Po-2{k,x )  = 


k  >  0, 


(21) 


G[k,x)  k\ 

where  for  each  k  >  1,  G(k,x)  is  the  c.d.f.  of  a  k-phase  Erlang  random  variable  with  parameter 
H  —  \q  and  G(0,  x)  =  1. 


Proof.  If  the  query  is  transmitted  to  the  (k  +  l)st  node,  then  it  had  k  successful  prior 
transmissions  without  expiring.  Let  T,;  denote  the  sojourn  time  at  the  ith  visited  node,  i  >  1. 
Because  each  node’s  transmission  queue  is  modeled  as  a  stable  M/M/1  queue,  {!)  :  i  >  1}  is  an 
i.i.d.  sequence  of  random  variables  with  parameter  p  —  Xq.  Denote  by  Y*.  the  total  time  elapsed  from 
the  moment  a  query  is  generated  at  an  uninformed  node  up  to  and  including  its  A;th  transmission, 
i.e., 

k 

U  =  £t., 

i= 1 

where  Y&  is  a  k-pliase  Erlang  random  variable  with  parameter  p  —  Xq.  It  is  well-known  (cf.  [23]) 
that,  for  k  >  1,  the  c.d.f.  of  is 

G{k,x)  =  P(Tfc  <  x)  =  1  -  V  ~  Xq^n  , 

t'o  ,,[ 

We  can  express  the  conditional  probability  Po2(k,x)  in  terms  of  the  random  variables  Y*.  and  Y^_ |_i 
by  noting  that 

P02(k,x )  =  P(Yjfc+i  >  x\Yk  <  x) 

is  the  probability  the  query  expires  at  the  kth  visited  node  while  awaiting  its  (k  +  l)st  transmission, 
given  it  had  successfully  made  k  prior  transmissions  and  was  active  just  before  joining  the  kth  node’s 
transmission  queue.  When  k  =  0,  Po2(0,x)  is  the  probability  the  query  expires  in  the  transmission 
queue  of  the  query  origin  node  given  by 

p02( 0,x)  =  P(Yi  >  X\X  =  x)=  P(Ti  >  x)  =  e~^~x^x. 
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For  k  >  1,  using  basic  conditional  probability, 


p02(k,x)  =  P(Yfc+i  >  x\Yk  <  x) 


G(k,x)  -G(k  +  l,x) 
G{k,x) 


e  [(/r  —  Ag)x]fc 

G(k,x)  k\ 


The  remaining  probabilities  in  row  0  of  P (k,x),  poo(k,x)  and  poi(k,x),  depend  on  whether  or  not 
the  query  revisits  uninformed  nodes  during  its  lifetime  when  r  <  oo.  For  this  reason,  it  is  necessary 
to  first  compute  the  probability  that  the  query  visits  a  particular  node  n  G  N  for  the  first  time  at 
its  kth  visit. 

To  this  end,  let  Uk  be  the  location  of  the  query  just  after  its  kth  hop  and  note  that  {Uk  :  k  >  0} 
is  a  time-homogeneous  DTMC  with  state  space  M  =  {1,...,JV}.  Define  its  one-step  transition 
probability  matrix  by  0(r )  =  As  in  section  2,  for  j  /  i,  let  p(i,j )  =  ||ajj  —  Xj  |  and  let 

di(r)  be  the  degree  of  node  i  G  AT.  Assuming  any  neighbor  of  the  current  node  is  equally  likely  to 
receive  a  query  transmission,  for  i,j  G  J\f  such  that  j  ^  the  transition  probability  9ij(r )  is 


1  /di(r),  if  p(i,j)<r, 
0,  if  p{i,j)>r. 


(Note  that  9u(r)  =  0  for  all  i  £  N  as  a  query  cannot  be  transmitted  to  the  current  node.) 

Now  let  q(k,  r )  be  the  probability  that  a  query  (or  event  agent)  visits  any  one  of  the  N  nodes 
for  the  first  time  at  the  kth  visit,  and  let  ur(i,j,k )  be  the  probability  of  visiting  node  j  at  least 
once  before  the  (k  +  l)st  visit,  given  that  the  query  (or  agent)  originates  at  node  i.  Let  wr(i,j,k) 
denote  the  probability  the  query  visits  state  j  for  the  first  time  on  the  kth  visit,  given  it  originated 
at  node  i.  We  have  the  following  important  lemma. 


Lemma  3  For  each  k  €  N  and  r  G  (0,  oo), 


q(k,  r)  ps  q(k,  r)  =  ^  ^  [ur(i,j,k)-ur(i,j,k-  1)] 


(22) 


where 


'i]\ 


(r)  +  ^2  dim  (r)ur(m,  j,  k  -  1),  k  >  1, 

m£Af\{j} 

0,  k  =  0. 


ur(i,j,k )  =  | 

Proof.  The  lemma  is  proved  using  standard  results  for  DTMCs.  Specifically,  define 


Tfj  =  inf {k  >  1  :  Uk  =  j\U0  =  i} 

as  the  first  passage  time  to  node  j  G  J\f,  given  that  the  query  (or  event  agent)  was  generated  at 
node  i  G  A f.  Then, 

ur(i,j,k)=P(T[j<k), 

and  these  probabilities  can  be  obtained  recursively  by  conditioning  on  the  location  of  the  query 
after  its  first  transmission.  The  derivation  is  similar  to  that  outlined  in  Theorem  4.1  of  [23]  and 
shows  that  for  k  >  1, 

ur(i,j,k)=9ij(r)+  ^2  0im(r)ur(m,j,  k  -  1),  i,jeAf, 
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where  ur(i,j,  0)  =  0  for  each  i.  j  £  J\f.  Using  ur(i,j,k),  the  probability  the  query’s  first  visit  to 
node  j  is  the  fcth  visit,  given  the  query  originated  at  node  i ,  is 


wr(i,j,k )  =  P (T[j  =  k)  =  ur(i,j,k )  -  ur(i,j,k  -  1),  k  >  1. 


Assuming  a  query  is  generated  at  any  i  £  N  with  equal  probability  (i.e.,  P(Uo  =  *)  =  1/-1V  for  all 
i  £  J\f),  via  unconditioning,  the  approximate  probability  a  query  visits  a  distinct  node  at  the  kth 
visit  is 


q(k, r )  ~  q(k, r ) 


JqY  Y  k>  1. 

i&M  j£Af\{i} 


Lemma  3  facilitates  simple  approximations  for  the  transition  probabilities  poo(k,x)  and  poi(k,x), 
k  >  0,  which  are  provided  in  the  next  proposition. 

Proposition  6  The  transition  probabilities  poo(k,x)  and  poi(k,x),  k>0,  are  respectively  approx¬ 
imated  by 

Poo(k,x )  «  [1  -  q(k  +  1,  r)(l  -  7r0(r))]  [1  -p02(k,x)\,  (23) 

Poi(k,  x)  «  q(k  +  l,r)[l  -  7 r0(r)]  [1  -p02(^,x)].  (24) 

Proof.  This  approximation  assumes  that  if  node  i  is  uninformed  when  a  query  first  visits 
the  node,  it  remains  uninformed  during  any  subsequent  visits  to  node  i  by  the  same  query.  We 
justify  this  assumption  by  noting  that  the  mean  recurrence  time  to  node  i  is  proportional  to  r.  To 
approximate  poo(k,x),  condition  on  whether  or  not  the  ( k  +  l)st  visited  node  is  distinct.  First, 
given  the  query  does  not  expire  at  the  kth.  visited  node,  the  ( k  +  l)st  visited  node  is  not  distinct 
with  probability  1  —  q(k  +  l,r).  In  the  second  case,  given  the  query  does  not  expire  at  the  kth 

visited  node,  the  ( k  +  l)st  node  is  distinct  with  probability  q{k  +  l,r),  and  it  is  uninformed  with 

probability  7To(r).  Therefore,  the  probability  of  locating  an  uninformed  node  at  the  (k  +  l)st  visit, 
given  the  query  was  active  just  before  joining  the  transmission  queue  of  fcth  node  is,  for  k  >  0, 

Poo{k,x)  «  [l-q(k  +  l,r)][l-po2(k,x)\  +  q(k  +  l,r)ir0(r)[l-p02(k,x)] 

=  [1  -  q(k  +  l,r)(l  -  7r0(r))]  [1  -po2(k,x)]. 

To  approximate  poi(k,x),  note  that  the  query  moves  from  state  0  (active)  to  state  1  (answered) 
if  it  was  successfully  transmitted  from  the  fcth  visited  node  to  a  distinct  node  that  is  informed. 
Therefore,  for  k  >  0, 

Poi{k,x)  «  q(k  +  l,r)[l  -  7r0(r)][l  -  p02(k,x)}. 


4.2  Range-Dependent  Query  Failure  Rate 

Using  the  approximation  of  P (k,x),  we  now  provide  improved  approximations  for  the  WSN 
traffic  rates,  the  steady  state  proportion  of  time  nodes  are  uninformed,  and  the  steady  state  pro¬ 
portion  of  failed  queries.  It  was  shown  in  Section  3  that,  if  each  sensor’s  range  is  such  that  all  N  —  1 
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other  nodes  belong  to  its  neighborhood,  the  total  arrival  rate  of  witnessed  events  (both  local  and 
external)  to  the  node’s  event  table  is 


A 


A  =  A 


'1  -  (1  -a)e+1' 
a 


The  approximation  A  does  not  account  for  the  revisiting  effects  noted  in  this  section.  The  following 
result  uses  q(k,r)  to  correct  for  revisits  and  improve  the  approximate  total  arrival  rate  to  the  event 
table.  To  distinguish  these  values,  let  A (r)  be  the  total  arrival  rate  of  local  and  external  events  as 
a  function  of  r.  Then  we  can  write 


A(r)  ~  A(r)  =  A  +  Ad(r) 


q(l,r)(l-a)  +  g(2,r)jl  -  a f  + _ ^q(£,r)(l  -  a)e 


d(r) 


d(r ) 


d{r ) 


=  A 


1  +  ^g(i,r)(! 


—  a 


i= 1 


where  d(r)  is  the  network’s  average  node  degree.  Using  A (r),  the  steady  state  proportion  of  time 
nodes  are  uninformed,  vro(r),  is 


no(r) 


exp 


A 

h 


l  +  ^2q(i,r)(l 


i= 1 


(25) 


Equation  (25)  is  used  to  compute  the  elements  of  P(fc,  x),  the  limiting  matrix  A(x),  and  ultimately 
the  limiting  probability  i>2(x)  to  obtain  A.r  via  (20).  The  asymptotic  validity  of  this  approximation 
is  discussed  in  the  next  subsection. 


4.3  Asymptotic  Validity  of  Approximation 

In  this  subsection,  we  show  that  the  finite  transmission  range  approximation  is  asymptotically 
valid  by  proving  that,  for  large  N ,  the  proportion  of  query  failures  converges  to  as  r  — >  oo.  To 
this  end,  we  have  the  following  important  lemma. 


Lemma  4  For  large  N,  as  r  — >•  oo,  q(k.  r)  — >  1  for  each  k  G  N. 
Proof.  First  note  that 


lim  dAr)  =  lim  N  1  (p(i,j)  <  r)  =  N  —  1. 

r—¥ oo  r—¥ oo  < 


Therefore,  for  i,j  €  A f  with  j  7^  i, 


Qij(r)  = 


di(r)  N  —  1 
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as  r  — v  oo.  By  induction  on  k  £  N,  we  now  characterize  the  limiting  behavior  of  ur(i,j,k)  as 
r  — V  oo.  For  k  =  1,  note  that  ur(i,j,  1)  =  9ij{r)  — V  1  /(N  —  1).  For  k  =  2,  it  is  easy  to  show  that 


lim  ur(i,j,  2) 

r— >•  oo 


lim  %(r)  +  V]  9im(r)ur(rn,j,  1) 

r— >-oo  \  z ' 

\  meA/'XO'} 

1  /  -i  \  2 

“TT +  E 


JV-  1 


m&Af\{i,j} 


1 


N  —  1 


N  -  1 


+  0(A 


-2\ 


where  0(iV  2)  — v  0  as  IV  — v  oo.  For  the  inductive  step,  assume  ur(i,j,n )  — V  n/(N  —  1)  +  0(N  2) 
for  any  n  £  N.  With  some  simplification  we  obtain 


lim  ur(i,j,n  + 1)  =  lim  %(r)  +  V'  9im(r)ur(m,j,n) 

r — Von  r — Von  l  '  < 

meAT\{j} 

r 

n 


i 


1+  E  v2: 


N  -1  ^  N  -  1 

m£Af\{i,j} 


N  —  1 


+  0(N 


-2-, 


n  +  1 
A-  1 


+  0(iY"2), 


which  completes  the  induction  proof.  Therefore,  for  each  k  £  N  and  i,j  £  Af  with  j  ^  i, 

1 


lim  wr(i,j,k )  =  lim  [ur(«,  j,k)  -  ur(i,j,k  -  1)]  = 

r— >-oo  r— >•  oo  iV  —  1 


and  consequently, 


lim  q(k,r)  =  lim  —  wr(i,j,k)  =  —  — — -  =  1. 

r- J-oo  V  ’  ’  r-KX)  iV  E^  E^  rv  ’ 1  N  E-/  E_^  AT  -  1 

iGA/”  j&A/"\{*}  i&A/-  jsA/”\{i} 


Lemma  4  will  now  be  used  to  prove  Theorem  2  which  asserts  that,  as  r  — V  oo,  the  approximate 
event  arrival  rate,  proportion  of  time  uninformed,  and  the  proportion  of  query  failures  all  converge 
appropriately  to  their  respective  infinite-range  counterparts  for  large  networks. 

Theorem  2  For  large  N,  as  r  — V  oo,  A(r)  —V  A,  7To(r)  —V  ttq,  and  Ar  — v  A^. 


Proof.  By  Lemma  4,  q(k,r)  —V  1  for  each  k  £  N  as  r  — V  oo.  Therefore, 


lim  A(r)  =  lim 

r— >oo  r—>  oo 


A  + Ay^g(z,r)(i 


—  a 


i= 1 


A  lim  )  q(i,r)(l 

r—±OG  E- — J 


;=o 


=  A 


1  -  (1  -  a 


y+ 1 


a 


—  a 


=  A. 


Consequently,  by  (25)  we  see  that  7To(r)  — V  7Tq  as  r  — V  oo.  Next,  recall  that  for  r  <  oo, 


A  ~  Ar  = 


V2(x)d  H(x) 
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where  ^(x)  =  lim*.^,*,  vk (x).  So  as  r  — )•  oo,  we  substitute  v°(x)  =  (tto,  1  —  vro,  0)  in  the  expression 

fc-i 

vk(x)  =  u°(x)  P (n,x). 


Using  (13),  (21),  (23),  and  (24),  we  now  show  by  induction  on  k  that  the  elements  of  vk(x)  are 

xo(x)  =  7rQ+1G(k,  x),  (26) 

k 

viix )  =  [^oC1  “  no)G{n,x)\  +  1  -7To,  (27) 

n=  1 

k- 1 

vUx)  =  1  -  7Tq_1G(A:,x)  -  (1  -  7T0)  ^7To_1G(n,x)  .  (28) 

71=1 

For  k  =  1,  applying  (14)  with  u°(x)  =  (7To,  1  —  7To,  0),  it  is  easy  to  see  that 

«o(®)  =  nlGil^x), 

v{(x)  =  7T0(1  -  7T0)G(l,x)  +  1  -  7T0, 

v\{x )  =  7T0[1  -  G(l,x)], 

where  the  summation  in  (28)  is  0  when  k  =  1.  Similarly,  for  k  =  2, 
vo(x)  = 

Ul(x)  =  7To(l  -  7T0)G(2,x)  +  7T0(1  -  7T0)G(l,x)  +  1  -  7T0, 
n|(x)  =  7T0  [1  -  (1  -  7T0)G(l,x)  -  7T0G(2,x)]  ; 

therefore,  the  result  holds  for  k  =  1,2.  For  the  inductive  step,  assume  that  (26)-(28)  hold  for  an 
arbitrary  m  £  N.  Then,  after  some  matrix  algebra,  we  obtain 

u™+1(x)  =  7r™+2G(tn  +  l,x), 

m+ 1 

V?+\x)  =  ^  [7To(l  -  7T0)G(n,x)]  +  1  -  7T0, 


(x)  =  7 r0  1  -  7rolG(?n  + l,x)  -  (1  -  7r0)^7TQ  1G(n,x)  , 


and  the  induction  proof  is  complete.  Now,  as  r  — >  oo, 


V2(x)  =  lim  vk{x)  =  lim  7To  (l  —  7Tq  1G(k,  x)  —  (1  —  7To)  tTq  1G(n,  x)^ 

k— too  k— too  \  J  I 

L  \  n=l  / J 


=  7T0 


1  -  (1  -7T0)^<  1G{n,x)  . 
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We  obtain  a  closed-form  expression  for  V2(x)  via  its  Laplace-Stieltjes  transform,  v2(s),  given  by 


v2 (s)  =  e  sxdv2(x )  =  7 r0 


!  _  1  -  Tip  /  7Tq(^ 

7 Tr\  ‘  \  //  —  ^ 


—  A„ 


=  7 T0 


=  7T0 


^0  “  V  M  -  N  +  s 

!  _  1  ~  7T0  /  7T0(/i  -  A, 

7 Trv  l  ^  \  //  — 


7T0 


^n=0 


\  —  Aq,  +  S 


-  1 


1  - 


(l-7r0)(/A-Ag) 

(1  -  7T0)(/i  -  Aq)  +  S_ 


Now,  V2  (s)  can  be  inverted  analytically  to  obtain 


v2(x)  =  /T 


-l  f  $2 (s) 


=  7 To  e-(l-«o)(A*-A9)*j 


where  £  1  is  the  inverse  Laplace  transform  operator.  Finally,  we  obtain 


lim  Ar  = 


t r0  e-(1-^o)^-A7)a:d//(x) 


P(T  >  A|/n  =  0,  A  =  x)7T0diL(x) 


=  P(T  >  X) 

=  Aoo. 


By  considering  a  limited  transmission  range,  the  model  presented  herein  explicitly  accounts 
for  query  revisiting  and  boundary  effects  (i.e. ,  the  effect  of  nodes  adjacent  to  the  boundary  of  the 
sensor  field).  In  Section  5,  we  illustrate  and  assess  the  quality  of  the  finite  and  infinite  transmission 
range  approximations  by  comparing  the  steady  state  proportion  of  time  uninformed  and  proportion 
of  query  failures  with  results  obtained  by  a  commercial  network  simulator. 

5  Numerical  Examples  and  Validation 

The  analytical  approximations  of  Sections  3  and  4  provide  a  relatively  easy  way  to  evaluate  the 
behavior  of  query-based  wireless  sensor  networks.  In  this  section,  we  assess  the  quality  and  validity 
of  these  approximations  by  comparing  them  with  simulated  values  obtained  using  the  OPNET 
commercial  network  simulator.  Presented  herein  are  summary  tables  and  figures  for  uniform- 
topology  networks  with  a  variety  of  distributional  assumptions  and  sensor  transmission  ranges.  For 
each  experiment,  the  minimum  transmission  range  was  chosen  to  ensure  a  connected  network  with 
probability  p  =  0.9999  using  (1).  Initially,  results  for  1000-  and  5000-node  networks  are  provided 
before  presenting  an  extensive  validation  study  that  examines  impact  of  our  model  assumptions. 

For  each  scenario,  the  performance  parameter  is  the  maximum  absolute  deviation  (MAD)  be¬ 
tween  the  approximated  value  and  its  simulated  counterpart  over  a  finite  set  of  time-to-live  values, 
L  =  {1,  2, . . . ,  30}.  We  choose  this  set  because,  for  many  typical  wireless  applications,  a  TTL 
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counter  between  3  and  25  is  suitable.  For  each  £  £  L,  let  7Tq  be  the  approximate  steady  state 
proportion  of  time  nodes  are  uninformed,  assuming  r  =  oo,  which  is  obtained  via  (6),  i.e., 


7r0  =  exp 


A  fl-(l-aY+1 
"6 


a 


Similarly,  let  7r(j(r)  be  the  same  value,  assuming  r  <  oo,  obtained  by  (25).  That  is, 

t  n 


7r0(r)  =  exp 


A 

"6 


1  +  X]5'(i,r)(1 


—  a 


i= 1 


For  both  cases,  the  probability  a.  is  approximated  using  the  fixed  point  algorithm  described  in 
Section  3.  To  express  the  dependence  of  A  on  the  TTL  value  £,  let  A^,  and  Af.  denote  the  steady 
proportion  of  query  failures  when  r  =  oo  and  r  <  oo,  respectively.  Using  (12)  and  (20),  respectively, 
we  compute 


Aqo  —  ^0 


and 


exp 


At  = 


-(1  -  TTq)(h  -  \q)x  d H(x) 


v2(x)d  H(x), 


where  v2(x)  is  obtained  via  (19).  In  cases  where  the  integrals  cannot  be  evaluated  in  closed  form, 
we  perform  numerical  integration  via  the  trapezoidal  rule.  Finally,  we  define  ttq(£)  as  the  simulated 
steady  state  proportion  of  time  nodes  are  uninformed,  and  As(£)  as  the  simulated  steady  state 
proportion  of  query  failures  when  the  TTL  counter  is  £  £  L. 

The  MAD  between  the  true  (simulated)  values  and  their  corresponding  analytical  approxima¬ 
tions  are  therefore 

Dn  =  max  |  7Tq  {£)  -  7?0  (£)  | ,  (29) 

where  7To(£)  =  ^  r  =  °°>  an^  tto(^)  =  if  r  <  oo.  Similarly,  let 


D/\  =  max 
t&L 


As{£)-Aq{£) 


(30) 


where  A(£)  =  A^  if  r  =  oo,  and  A{£)  =  Af  if  r  <  oo.  For  Examples  1  and  2  that  follow,  a 
few  parameter  values  were  held  constant;  these  values  are  summarized  in  Table  1.  Moreover,  we 
assumed  event  lifetimes  are  exponentially  distributed  with  parameter  5  in  these  two  cases,  but  this 
assumption  is  relaxed  in  Example  3. 


Table  1:  Summary  of  parameter  values  for  OPNET  simulation:  Examples  1  and  2. 


Parameter 

Parameter  description 

Value 

Transmitter’s  exponential  transmission  rate 

5.000 

A 

Poisson  rate  of  locally-witnessed  events  (for  all  n  £  M) 

0.005 

7 

Poisson  rate  of  locally-generated  queries  (for  all  n  £  Af) 

0.050 

1/5 

Mean  event  lifetime 

10.000 

1/(5 

Mean  query  lifetime  (for  all  distributions) 

5.000 
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The  analytical  approximations  were  coded  in  the  C  programming  language  and  executed  in 
Microsoft®  Visual  Studio®  2008  on  a  personal  computer  equipped  with  an  Intel®  Core™  2  Duo 
CPU  operating  at  3.00GHz  with  2.00  GB  of  RAM.  The  simulated  values  were  obtained  via  a 
discrete-event  simulation  model  created  in  the  OPNET  Modeler  Wireless  Suite  v.  15.  Ten  (10) 
independent  replications  were  performed  for  each  i  €  L  to  ensure  a  standard  error  less  than  5  x  10~4. 
The  reported  simulated  values  represent  the  average  of  the  10  replications.  The  run  length  was 
3720s,  including  a  120s  warm-up  period  for  each  replication.  The  simulation  experiments  were 
conducted  on  a  personal  computer  equipped  with  an  Intel®  Core™  i7  CPU  operating  at  2.67GHz 
with  2.00  GB  of  RAM. 

Example  1:  1000-Node  Network:  Here,  we  present  results  for  a  1000-node  wireless  sensor 
network  with  nodes  distributed  randomly  in  a  3335m  x  3335m  sensor  field.  The  node  density  is 
'(/’  ~  9.00  x  10-5  nodes  per  square  meter.  To  ensure  a  connected  network  with  probability  0.9999,  the 
minimum  required  sensor  transmission  range  is  r  =  239m.  Therefore,  we  considered  the  following 
transmission  ranges:  350m,  500m,  1000m,  5000m.  By  doing  so,  we  are  able  to  assess  the  quality  of 
both  the  infinite  and  finite  range  approximations.  Table  2  summarizes  the  MAD  in  the  proportion 
of  time  uninformed  using  each  transmission  range.  The  column  labeled  “r  =  oo”  corresponds  to 
the  infinite  transmission  range  approximation,  and  the  column  labeled  “r  <  oo”  is  the  finite  range 
approximation. 


Table  2:  MAD  in  the  proportion  of  time  uninformed  ( D when  N  =  1000. 


Query  lifetime 

350m 

r  =  oo  r  <  oo 

500m 

r  =  oo  r  <  oo 

1000m 

r  =  oo  r  <  oo 

5000m 

r  =  oo  r  <  oo 

Exponential(0.2) 

0.0523 

0.0045 

0.0289 

0.0065 

0.0202 

0.0116 

0.0024 

0.0060 

Triangular(0.1,  5.0,  9.9) 

0.0529 

0.0059 

0.0296 

0.0055 

0.0173 

0.0088 

0.0173 

0.0131 

Uniform(0.1,  9.9) 

0.0507 

0.0068 

0.0285 

0.0041 

0.0171 

0.0093 

0.0039 

0.0070 

Rayleigh(5.645) 

0.0511 

0.0050 

0.0298 

0.0055 

0.0176 

0.0088 

0.0043 

0.0069 

Weibull(3.0,  5.6) 

0.0511 

0.0061 

0.0302 

0.0061 

0.0173 

0.0088 

0.0035 

0.0072 

Table  2  indicates  an  order  of  magnitude  improvement  by  including  the  revisiting  effect,  especially 
when  the  actual  transmission  range  is  small  (350m).  Because  queries  are  more  likely  to  revisit 
neighbors  when  the  transmission  range  is  small,  the  difference  between  the  two  approximations  is 
quite  pronounced.  Table  2  also  illustrates  the  performance  of  the  approximations  when  the  query 
lifetime  distribution  is  not  nrenroryless.  The  results  are  consistent  with  the  case  of  exponential 
query  lifetimes.  Figure  3  shows  the  performance  of  the  approximations  and  reveals  that  the  finite 
range  approximation  is  superior  to  the  infinite  range  approximation  for  all  TTL  values  when  r  is 
small.  Indeed,  the  gap  between  the  latter  approximation  and  OPNET  simulation  values  increases 
with  £  since  the  revisiting  effect  is  more  pronounced  when  the  time-to-live  counter  is  large.  For 
larger  ranges,  the  approximations  nearly  coincide  and  both  closely  track  the  simulated  values. 

Results  for  the  steady  state  proportion  of  query  failures  are  summarized  in  Table  3.  Both 
approximation  schemes  perform  extremely  well  (the  maximum  absolute  deviation  over  all  cases 
is  less  than  0.049).  It  is  also  worth  noting  that  the  finite  range  approximation  outperforms  the 
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350m  Transmission  Range  500m  Transmission  Range 


1 000m  T  ransmission  Range  5000m  T  ransmission  Range 


Figure  3:  Comparison  of  7ro  values  with  Weibull  query  lifetimes  (N  =  1000):  (-)  OPNET;  (o)  r  =  oo;  (+)  r  <  oo. 

infinite  range  approximation,  particularly  when  r  is  relatively  small.  The  results  are  consistent  for 
non-memoryless  query  lifetimes.  Figure  4  graphically  depicts  the  four  cases. 


Table  3:  MAD  in  the  proportion  of  failed  queries  (Da)  when  N  =  1000. 


Query  lifetime 

350m 

r  =  oo  r  <  oo 

500m 

r  =  oo  r  <  oo 

1000m 

r  =  oo  r  <  oc 

5000m 

r  =  oo  r  <  oo 

Exponential(0.2) 

0.0371 

0.0246 

0.0168 

0.0103 

0.0047 

0.0055 

0.0060 

0.0052 

Triangular(0.1,  5.0,  9.9) 

0.0459 

0.0301 

0.0170 

0.0128 

0.0030 

0.0008 

0.0082 

0.0024 

Uniform(0.1,  9.9) 

0.0383 

0.0258 

0.0178 

0.0121 

0.0035 

0.0015 

0.0051 

0.0016 

Rayleigh(5.645) 

0.0237 

0.0127 

0.0065 

0.0158 

0.0256 

0.0270 

0.0255 

0.0247 

Weibull(3.0,  5.6) 

0.0485 

0.0306 

0.0230 

0.0149 

0.0031 

0.0014 

0.0022 

0.0014 
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350m  Transmission  Range 


Time  to  Live  Counter 


Time  to  Live  Counter 


5000m  Transmission  Range 


Figure  4:  Comparison  of  A  values  with  Weibull  query  lifetimes  (N  =  1000):  (-)  OPNET;  (o)  r  =  oo;  (+)  r  <  oo. 


Example  2:  5000-Node  Network:  Here,  we  consider  a  5000-node  wireless  sensor  network  with 
nodes  deployed  in  the  same  region  as  the  1000- node  case  but  with  node  density  ~  4.50  x  10~4  nodes 
per  square  meter.  To  ensure  a  connected  network  with  probability  0.9999,  the  minimum  required 
sensor  transmission  range  is  r  =  112m.  Therefore,  we  considered  the  following  transmission  ranges: 
115m,  350nr,  500m,  and  5000m.  Table  4  illustrates  the  quality  of  both  approximations  for  the 
5000-node  network.  The  maximum  absolute  deviation  for  the  proportion  of  time  uninformed  is  less 
than  0.082  for  r  =  oo,  and  it  is  reduced  to,  at  most,  0.0174  when  the  revisiting  effect  is  included.  As 
before,  the  superiority  of  the  finite  range  approximation  is  generally  more  pronounced  for  smaller 
transmission  ranges.  Figure  5  depicts  the  simulated  and  approximated  values  of  7Tq  when  the  query 


Table  4:  MAD  in  the  proportion  of  time  uninformed  (Dv)  when  N  =  5000. 


Query  lifetime 

115m 

r  =  oo  r  <  oo 

350m 

r  =  oo  r  <  oo 

500m 

r  =  oo  r  <  oo 

5000m 

r  =  oo  r  <  oo 

Exponential(0.2) 

0.0605 

0.0172 

0.0107 

0.0019 

0.0082 

0.0031 

0.0040 

0.0049 

Triangular(0.1,  5.0,  9.9) 

0.0612 

0.0174 

0.0100 

0.0015 

0.0068 

0.0026 

0.0053 

0.0062 

Uniform(0.1,  9.9) 

0.0819 

0.0054 

0.0105 

0.0013 

0.0074 

0.0025 

0.0048 

0.0058 

Rayleigh(5.645) 

0.0595 

0.0172 

0.0101 

0.0016 

0.0074 

0.0026 

0.0051 

0.0060 

Weibull(3.0,  5.6) 

0.0611 

0.0174 

0.0099 

0.0014 

0.0068 

0.0028 

0.0073 

0.0083 
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lifetime  follows  a  triangular  distribution.  When  the  transmission  range  is  small  (115m),  we  see 
some  discrepancy  between  the  two  approximation  schemes.  However,  for  the  other  three  cases,  the 
approximations  nearly  coincide  and  are  very  similar  to  the  simulated  results  ( Dn  <  0.011). 

1 1 5m  T ransmission  Range  350m  T ransmission  Range 


500m  Transmission  Range  5000m  Transmission  Range 


Figure  5:  Comparison  of  no  values  with  triangular  query  lifetimes  (N  =  5000):  (-)  OPNET;  (o)  r  =  oo;  (+)  r  <  oo. 


Next,  we  compare  the  maximum  absolute  deviation  of  the  proportion  of  query  failures.  Table 
5  shows  that  the  maximum  deviation  values  are  bounded  above  by  0.0725.  Again,  the  finite 
range  approximation  outperforms  the  infinite  range  version  when  the  transmission  range  is  small. 
However,  for  larger  ranges,  the  results  nearly  coincide  and  closely  track  the  simulated  values. 


Table  5:  MAD  in  the  proportion  of  query  failures  (D a)  when  N  =  5000. 


Query  lifetime 

115m 

r  =  oo  r  <  oo 

350m 

r  =  oo  r  <  oo 

500m 

r  =  oo  r  <  oo 

5000m 

r  =  oo  r  <  oo 

Exponential(0.2) 

0.0493 

0.0283 

0.0046 

0.0022 

0.0045 

0.0043 

0.0061 

0.0044 

Triangular(0.1,  5.0,  9.9) 

0.0588 

0.0333 

0.0052 

0.0026 

0.0013 

0.0024 

0.0044 

0.0045 

Uniform(0.1,  9.9) 

0.0724 

0.0493 

0.0044 

0.0047 

0.0051 

0.0021 

0.0078 

0.0023 

Rayleigh(5.645) 

0.0371 

0.0170 

0.0214 

0.0192 

0.0265 

0.0225 

0.0286 

0.0228 

Weibull(3.0,  5.6) 

0.0619 

0.0351 

0.0071 

0.0064 

0.0025 

0.0041 

0.0039 

0.0021 

Figure  6  graphically  depicts  the  simulated  and  approximated  values  of  A  and  illustrates  the 
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high  quality  of  the  approximations.  In  the  worst  case  (115m),  the  MAD  is  less  than  0.0725  and 
0.05  for  r  =  oo  and  r  <  oo,  respectively. 


115m  Transmission  Range 


350m  Transmission  Range 


500m  Transmission  Range  5000m  Transmission  Range 


Figure  6:  Comparison  of  A  values  with  triangular  query  lifetimes  ( N  =  5000):  (-)  OPNET;  (o)  r  =  oo;  (+)  r  <  oo. 


Example  3:  Model  Validation:  Finally,  we  conducted  an  experiment  to  validate  the  approxima¬ 
tions  when  some  of  the  model  assumptions  are  violated.  For  the  benchmark  simulation  experiments 
presented  here,  events  arrive  according  to  a  renewal  process  with  a  specified  (non-exponential)  in¬ 
terarrival  time  distribution  (i.e.,  the  event  arrival  process  is  not  Poisson).  This  experiment  also 
employs  non-exponential  event  agent  and  query  lifetimes,  both  of  which  are  incorporated  in  the 
approximations  of  Sections  3  and  4. 

Table  6  provides  a  summary  of  the  numerical  results  for  45  distinct  test  cases  using  a  1000-node 
wireless  sensor  network  with  nodes  distributed  randomly  in  a  3335m  x  3335m  sensor  field.  The 
node  density  is  ip  ~  9.00  x  10~5  nodes  per  square  meter.  To  ensure  a  connected  network  with 
probability  0.9999,  the  minimum  required  sensor  transmission  range  is  r  =  239m;  therefore,  we  set 
r  =  350m. 
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Table  6:  MAD  in  the  proportion  of  time  uninformed  (Dv)  and  proportion  of  failed  queries  (Da)- 


Trial 

Event  interarrival  time 

Event  lifetime 

Query  lifetime 

A 

r  =  oo 

T 

r  <  oo 

Da 

r  =  oo  r  <  oo 

1 

Erlang(5,  40.0) 

Erlang(4,  2.5) 

Erlang(5,  1.0) 

0.0711 

0.0233 

0.0421 

0.0265 

2 

Erlang(5,  40.0) 

Triangular(0.1.  10.0, 

19.9) 

Erlang(5,  1.0) 

0.0693 

0.0215 

0.0428 

0.0273 

3 

Erlang(5,  40.0) 

Uniform(0.1.  19.9) 

Erlang(5,  1.0) 

0.0596 

0.0124 

0.0424 

0.0273 

4 

Erlang(5,  40.0) 

Erlang(4,  2.5) 

Rayleigh(5.645) 

0.0696 

0.0213 

0.0201 

0.0106 

5 

Erlang(5,  40.0) 

Triangular(0.1,  10.0, 

19.9) 

Rayleigh(5.645) 

0.0699 

0.0220 

0.0203 

0.0106 

6 

Erlang(5,  40.0) 

Uniform(0.1.  19.9) 

Rayleigh(5.645) 

0.0594 

0.0130 

0.0213 

0.0106 

7 

Erlang(5,  40.0) 

Erlang(4,  2.5) 

Triangular(0.1,  5.0,  9.9) 

0.0687 

0.0213 

0.0399 

0.0262 

8 

Erlang(5,  40.0) 

Triangular(0.1.  10.0, 

19.9) 

Triangular(0.1,  5.0,  9.9) 

0.0682 

0.0206 

0.0407 

0.0267 

9 

Erlang(5,  40.0) 

Uniform(0.1.  19.9) 

Triangular(0.1,  5.0,  9.9) 

0.0584 

0.0100 

0.0400 

0.0259 

10 

Erlang(5,  40.0) 

Erlang(4,  2.5) 

Uniform (0.1,  9.9) 

0.0691 

0.0218 

0.0340 

0.0209 

11 

Erlang(5,  40.0) 

Triangular(0.1.  10.0, 

19.9) 

Uniform (0.1,  9.9) 

0.0690 

0.0211 

0.0342 

0.0205 

12 

Erlang(5,  40.0) 

Uniform(0.1.  19.9) 

Uniform (0.1,  9.9) 

0.0589 

0.0118 

0.0344 

0.0211 

13 

Erlang(5,  40.0) 

Erlang(4,  2.5) 

Weibull(3.0,  5.6) 

0.0699 

0.0229 

0.0424 

0.0278 

14 

Erlang(5,  40.0) 

Triangular(0.1.  10.0, 

19.9) 

Weibull(3.0,  5.6) 

0.0677 

0.0201 

0.0427 

0.0282 

15 

Erlang(5,  40.0) 

Uniform(0.1.  19.9) 

Weibull(3.0,  5.6) 

0.0584 

0.0096 

0.0428 

0.0278 

16 

Triangular (1.0,  200.0, 

399.0) 

Erlang(4,  2.5) 

Erlang(5,  1.0) 

0.0739 

0.0261 

0.0438 

0.0282 

17 

Triangular(1.0,  200.0, 

399.0) 

Triangular(0.1.  10.0, 

19.9) 

Erlang(5,  1.0) 

0.0713 

0.0235 

0.0429 

0.0274 

18 

Triangular (1.0,  200.0, 

399.0) 

Uniform(0.1.  19.9) 

Erlang(5,  1.0) 

0.0615 

0.0145 

0.0433 

0.0282 

19 

Triangular (1.0,  200.0, 

399.0) 

Erlang(4,  2.5) 

Rayleigh(5.645) 

0.0718 

0.0235 

0.0213 

0.0106 

20 

Triangular (1.0,  200.0, 

399.0) 

Triangular(0.1.  10.0, 

19.9) 

Rayleigh(5.645) 

0.0728 

0.0250 

0.0215 

0.0106 

21 

Triangular(1.0,  200.0, 

399.0) 

Uniform(0.1.  19.9) 

Rayleigh(5.645) 

0.0614 

0.0150 

0.0217 

0.0114 

22 

Triangular(1.0,  200.0, 

399.0) 

Erlang(4,  2.5) 

Triangular(0.1,  5.0,  9.9) 

0.0713 

0.0239 

0.0410 

0.0273 

23 

Triangular(1.0,  200.0, 

399.0) 

Triangular(0.1.  10.0, 

19.9) 

Triangular(0.1,  5.0,  9.9) 

0.0725 

0.0249 

0.0405 

0.0265 

24 

Triangular(1.0,  200.0, 

399.0) 

Uniform(0.1.  19.9) 

Triangular(0.1,  5.0,  9.9) 

0.0523 

0.0052 

0.0418 

0.0277 

25 

Triangular(1.0,  200.0, 

399.0) 

Erlang(4,  2.5) 

Uniform (0.1,  9.9) 

0.0709 

0.0237 

0.0343 

0.0212 

26 

Triangular(1.0,  200.0, 

399.0) 

Triangular(0.1.  10.0, 

19.9) 

Uniform (0.1,  9.9) 

0.0717 

0.0238 

0.0356 

0.0219 

27 

Triangular (1.0,  200.0, 

399.0) 

Uniform(0.1.  19.9) 

Uniform (0.1,  9.9) 

0.0618 

0.0146 

0.0296 

0.0163 

28 

Triangular (1.0,  200.0, 

399.0) 

Erlang(4,  2.5) 

Weibull(3.0,  5.6) 

0.0718 

0.0247 

0.0437 

0.0290 

29 

Triangular (1.0,  200.0, 

399.0) 

Triangular(0.1.  10.0, 

19.9) 

Weibull(3.0,  5.6) 

0.0721 

0.0246 

0.0431 

0.0286 

30 

Triangular (1.0,  200.0, 

399.0) 

Uniform(0.1.  19.9) 

Weibull(3.0,  5.6) 

0.0519 

0.0037 

0.0438 

0.0288 

31 

Uniform(1.0,  399.0) 

Erlang(4,  2.5) 

Erlang(5,  1.0) 

0.0737 

0.0259 

0.0453 

0.0297 

32 

Uniform(1.0,  399.0) 

Triangular(0.1.  10.0, 

19.9) 

Erlang(5,  1.0) 

0.0675 

0.0201 

0.0459 

0.0304 

33 

Uniform(1.0,  399.0) 

Uniform(0.1.  19.9) 

Erlang(5, 1.0) 

0.0630 

0.0158 

0.0457 

0.0306 

34 

Uniform(1.0,  399.0) 

Erlang(4,  2.5) 

Rayleigh(5.645) 

0.0742 

0.0260 

0.0243 

0.0127 

35 

Uniform (1.0,  399.0) 

Triangular(0.1.  10.0, 

19.9) 

Rayleigh(5.645) 

0.0737 

0.0258 

0.0234 

0.0123 

36 

Uniform (1.0,  399.0) 

Uniform(0.1.  19.9) 

Rayleigh(5.645) 

0.0637 

0.0167 

0.0233 

0.0133 

37 

Uniform (1.0,  399.0) 

Erlang(4,  2.5) 

Triangular(0.1,  5.0,  9.9) 

0.0751 

0.0276 

0.0435 

0.0298 

38 

Uniform (1.0,  399.0) 

Triangular(0.1.  10.0, 

19.9) 

Triangular(0.1,  5.0,  9.9) 

0.0745 

0.0269 

0.0440 

0.0300 

39 

Uniform (1.0,  399.0) 

Uniform(0.1.  19.9) 

Triangular(0.1,  5.0,  9.9) 

0.0634 

0.0155 

0.0444 

0.0300 

40 

Uniform (1.0,  399.0) 

Erlang(4,  2.5) 

Uniform (0.1,  9.9) 

0.0749 

0.0277 

0.0369 

0.0255 

41 

Uniform (1.0,  399.0) 

Triangular(0.1.  10.0, 

19.9) 

Uniform(0.1,  9.9) 

0.0749 

0.0270 

0.0356 

0.0219 

42 

Uniform (1.0,  399.0) 

Uniform(0.1.  19.9) 

Uniform (0.1,  9.9) 

0.0634 

0.0165 

0.0379 

0.0252 

43 

Uniform (1.0,  399.0) 

Erlang(4,  2.5) 

Weibull(3.0,  5.6) 

0.0728 

0.0257 

0.0458 

0.0311 

44 

Uniform (1.0,  399.0) 

Triangular(0.1.  10.0, 

19.9) 

Weibull(3.0,  5.6) 

0.0746 

0.0270 

0.0464 

0.0319 

45 

Uniform(1.0,  399.0) 

Uniform(0.1.  19.9) 

Weibull(3.0,  5.6) 

0.0648 

0.0161 

0.0469 

0.0319 

Table  6  reveals  some  very  interesting  results.  First,  we  note  that  the  performance  of  the  finite-range 
approximation  is  similar  to  that  reported  in  Example  1  which  assumed  Poisson-generated  events. 
Despite  the  fact  that  the  event  arrival  process  is  distinctly  non-Poisson,  and  the  query  and  event 
lifetimes  are  non-memory  less,  the  proportion  of  failed  queries  is  approximated  very  closely  using  the 
finite-range  approximation  (as  compared  to  the  OPNET  simulated  values).  Over  all  45  test  cases, 
the  maximum  absolute  deviation  is  on  the  order  of  0.03.  These  findings  are  significant  because 
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they  provide  empirical  evidence  that  the  approximations  are  not  heavily  influenced  by  the  Poisson 
arrival  assumption  at  the  event  tables  or  the  transmission  queues.  Moreover,  for  each  instance,  the 
approximated  proportion  of  time  uninformed  and  proportion  of  query  failures  ions  were  computed 
in  less  than  20  minutes  as  opposed  to  the  OPNET  simulation  results,  which  required  approximately 
2  hours. 
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